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ABSTRACT 


uroup  6~  has  designed  and  built  an  isoplanometer.  an  instrument  which  remotely  senses  the  isoplanatic  angle 
associated  with  short-term  atmospheric  conditions.  Verification  of  the  instrument  required  nighttime  testing 
using  stellar  sources  to  measure  current  atmospheric  condition.-*.  Preliminary  results  from  a  verification 
experiment  and  a  brief  description  of  the  instrument  are  presented.  Full  instrument  verification  requires 
locating  the  instrument  at  a  better  astronomical  site:  therefore,  our  data  should  not  be  considered  a  valid 
characterization  of  the  atmospheric  conditions  in  Lexington.  Massachusetts. 
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PRELIMINARY  EXPERIMENTAL  MEASUREMENT 
OF  ISOPLANATIC  ANGLE 


I.  INTRODUCTION 

Group  67  has  designed  and  built  an  isoplanometer.  an  instrument  which  remotely  measures  the 
isoplanatic  angle  associated  with  atmospheric  turbulence  conditions.  The  isoplanatic  angle  is  loose!)  defined 
as  the  angular  field  over  which  the  spatial  transfer  function  of  the  atmosphere  can  be  assumed  to  be  invariant. 
The  theory  and  design  of  the  isoplanometer  are  detailed  ir.  a  separate  document. 1  Here,  a  brief  rev  iew  ot  the 
theory  and  a  description  of  the  instrument  constructed  are  given.  Isoplanatic  angle  data  are  presented  for  a 
single  measurement  session,  taken  from  the  A-building  roof  of  Lincoln  Laboratory  on  25  July.  19XN. 
Measurement  procedures  and  instrument  verification  tests  also  are  discussed.  This  report  documents  the 
operation  of  the  instrument  (which  was  designed  to  operate  only  under  "good"  atmospheric  seeing 
conditions),  but  does  not  document  atmospheric  conditions  in  Lexington.  Massachusetts.  These  preliminary 
measurement^  should  not  be  used  as  verified  data. 

II.  THEORY 

For  completeness,  we  rev  iew  the  theory  used  for  these  measurements.  The  basis  of  the  theory  is  the 
relationship  ( first  suggested  by  Loos  and  Hogge~ )  betw  een  the  observed  intensity  fluctuations  of  starlight  and 
the  isoplanatic  angL.  The  resulting  measurement  technique  is  indirect. '  a  crucial  point  since  a  direct  mea¬ 
surement  would  require  a  large  amount  of  airborne  equipment  at  considerably  greater  complication  and 
expense. 

Consider  a  telescope  whose  aperture  has  been  apodized  by  a  circularly  symmetric,  spatially  varying 
intensity  mask.  Let  the  intensity  transmittance  of  this  mask  be  defined  by  P(p),  where  0  <  P(p)  <  1 .  and  p  is 
the  radial  distance  to  the  center  of  the  mask.  The  optical  power  collected  from  a  star  viewed  through  this 
apodized  telescope  is  denoted  as  S.  Atmospheric  turbulence  causes  starlight  intensity  to  vary  (twinkle): 
therefore.  S  is  time  variant.  It  can  be  show  n  that  the  normalized  variance  of  the  received  intensity  fluctuations 
of  a  star  observed  through  ths  apodized  telescope  aperture  is  given  by1 

o;iO/./j  =  4(  27n4().033secK'  '(07)k2  A""  |  CpizlWuidz  tli 


where 

A  =  2n  !  p  Pi p  i  dp 

k  =  2k/K  (A  =  viewing  wavelength i 

6  =  viewed  zenith  ancle  of  the  star 
/ 

z  =  slant  range  distance  to  the  star 

C"(z)  =  refractive  structure  parameter  of  the  atmosphere  measured  along  the  slant  range 
path  to  the  star 


I 


oc  7  r  1 2 

W(z)  =  J“|  JdppJ0(Lp)P(Lp)j  L“8/-’sin2  ~  dL 

Equation  ( 1 1  explicit!)  indicates  that  the  normalized  variance  is  a  function  of  both  zenith  viewing  angle 
<i>  ,  and  viewing  waselength  X.  (k  =  2tt/X ).  The  wavelength  dependence,  however,  is  a  small  second-order 
effect,  and  thus  the  wavelength  dependence  of  cr“(4>  )  will  be  dropped. 

Equation  ( 1 )  was  derived  under  the  assumption  of  weak  scattering,  but  the  result  is  known  to  remain  valid 
even  outside  this  regime.2'6  There  is  some  cr*  value,  however,  beyond  w  hich  Equation  ( 1 )  is  no  longer  correct. 
This  maximum  value  is  a  function  of  the  mask  P(p),  and  its  determination  is  an  area  of  active  research.  When 
ct~  is  less  than  the  allow  ed  maximum,  the  normalized  variance  is  said  to  be  unsaturated,  and  the  results  of  the 

S 

isoplanatic  angle  approximation  presented  below  are  considered  valid. 

T^e  isoplanatic  angle  B>  measured  along  the  slant  path  z  is  given  by"*' 

80(<5?.A  =  0.527  k~6/5 

The  zenith-angle  dependence  is  inherent  in  the  slant  path  z.  If  an  apodizing  mask  P(p)  can  be  chosen  such 
that  W(z)  =  cz5  \  then  it  follows  from  Equations  ( 1 )  and  (2)  that  the  observed  starlight  intensity  fluctuations 
can  be  related  to  the  isoplanatic  angle  by 

6o(0z.>s)  =  12.9 sec-8/5(0z)  A”6/5c?/5[o2(0)]  ‘V?  (3, 

which  is  the  desired  result. 

A  mask  satisfying  W(z|  =  cz5/'  at  a  wavelength  cf  X  =  5000  A  was  found  by  an  exhaustive  search 
algorithm.  The  search  was  restricted  to  binary  masks  [i.e..  P(  p)  =  0  or  1],  since  they  are  easier  to  fabricate. 
The  masks  examined  were  required  to  fit  on  an  8-in-diam.  Celestron  Schmidt-Cassegrain  telescope  to  be  user* 
in  the  experiment.  The  telescope  has  a  2. 8-in-diam.  secondary' mirror  obscuration;  therefore.  P(p)  =  0for 
p  <  1 .4  in.  Figure  1  shows  the  mask  solution  obtained  by  exhaustive  search.  The  black  areas  of  the  mask  are 
totally  obscuring:  the  white  areas  are  assumed  to  be  100-percent  transmissive.  The  radii  of  the  optimum  mask 
are  as  follow  s: 

Central  obscuration  1.4-in  radius 

Telescope  diameter  4.0  in  radius 

Innermost  clear  aperture  From  1.472-  to  1.726-in  radii 

Second  clear  aperture  From  2.476-  to  2.726-in  radii 

Outermost  clear  aperture  From  3.226-  to  4. 000-in  radii 


rv 
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(z)z5/3dz 


radians 
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Fit; wc  i  Mask  design 

Figure  2  shows  the  corresponding  percentage  error  in  approximating  the  7/  '  profile  b>  ad)usting  Wi/i. 
plotting  [Wi7)/e  *  z'  '-!]*  100.  The  value  of  A  for  this  mask  is  0.0150  nr.  and  the  value  for  c  is  X.S4"  • 

1  O' 1  nT.  As  seen  in  Figure  2.  the  approximation  error  is  quite  small.  At  zenith  angles  less  than  60".  and  under 
all  but  the  most  sev  ere  atmospheric  conditions,  contributions  to  Equation  (2)  are  insignificant  outside  0.5  to 
40  km.  Errors  are  less  than  +3  percent  over  the  most  crucial  0.5-  to  20-km  range,  where  the  C~  profile  is  most 
significant  under  “good"  atmospheric  conditions. 

Plugging  in  constants,  for  this  mask  the  relation  between  isoplanalic  angle  and  starlight  intensity  variance 
is 

^  5  r  1—3/  ^ 

60i vrk)  =  0. 443 1 46 sec  M,'(O7)  02(O)  pradians  (4i 

Note  that  the  zenith-angle  dependence  of  the  isoplanatic  angle  here  assumes  that  the  C ‘  profile  of  the 
atmosphere  is  on!  v  a  function  of  height,  or  that  conditions  are  locally  constant.  For  poor  atmospheric  sites 
(such  as  Lexington),  this  assumption  can  be  invalid. 

Walters  has  developed  a  technique  for  verifying  proper  isoplanometer  operation. 4  Data  are  collected  bv 
rapidiy  (even  few  minutes  1  changing  between  stars  w  ith  different  zenith  angles.  Equation  ( 1  )  indicates  that 
a  log-log  plot  of  normalized  variance  v  s  zenith  angle,  as  shown  by  Walters,  should  be  linear  u  ith  a  slope  of 
S/3.  The  degree  to  w  hich  the  measured  data  points  fit  this  linear  plot  sen  es  as  a  test  of  proper  isoplanometer 
operation.  This  is  valid,  howev  er.  onlv  if  ( 1 1  the  normalized  variance  is  not  saturated,  and  ( 2  >  w  hen  C  ~  is  onlv 

n 

a  function  of  height.  It  the  atmospheric  conditions  are  constant  over  long  periods,  an  equivalent  method  is 
to  track  the  same  star  as  zenith  angles  change:  however,  “constant"  atmosphere  is  rarelv  the  case,  even  for 
good  astronomical  site^ 


HEIGHT  (km) 
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III.  THK  INS  I  HI  MI  N  I 

The  experimental  setup  is  shown  in  Figure  3.  The  primary  components  consist  o|  a  C'elestron  s-m  f  I  " 
Schmidt-Cassegrainian  telescope,  an  apodi/.mg  mask,  an  optical  detector  assembly.  a  high-voltage  powc: 
supply.  processing  electronics,  an  oscilloscope,  and  a  computer  controller.  The  optical  detector  assembly 
consists  of  a  black  anodized  aluminum  tube  which  contains  a  beamsplitter,  a  25-mm  ret ic led  eyepiece.  a 
narrowband  filter,  and  a  photomultiplier  tube,  as  shown  in  Figure  4.  Starlight  passes  through  the  niasi 
through  the  telescope,  and  hits  the  beamsplitter.  A  portion  of  the  light  is  split  off  to  the  eyepiece,  which  o 
used  for  alignment  and  tracking  purposes.  The  exact  amount  of  light  split  off  is  not  important,  since  the 
isoplanatie  angle  calculation  is  independent  of  the  absolute  signal  level.  The  remaining  light  passes  through 
a  100-A-wide  optical  filter  centered  at  5000  A.  The  resulting  light  is  placed  on  a  photomultiplier  tube  t  P\FI  i 
In  order  to  increase  the  image  size  on  the  tube,  the  PMT  is  placed  nominally  at  1  in  behind  the  back  local  plane 
of  the  telescope.  Lsing  a  larger  than  diffraction-limited  spot  reduces  the  effects  of  PMT  nonuniformities.  The 
PMT  (an  Oriel  77340)  is  powered  by  an  Oriel  model  7070  high-voltage  power  supplv  and  current  momtoi 

the  electronics  unit  was  designed  by  E.  Coiagiuri  ;c  compute  certain  statistics  of  the  output  current  signal 
from  the  PMT.  It  consists  of  a  transimpedance  amplifier,  followed  by  a  voltage  amplifier,  variable-bandw  ldth 
5-pole  low -pass  filter,  12-bit  A/D  converter,  and  interlace  logic.  The  A/D  samples  the  resulting  signal  ai  a 
1-kHz  rate,  and  digital  logic  computes  the  sum  of  the  samples  and  the  sum  of  the  square  of  the  samples.  The 
bandwidth  of  the  scintillations  rolled  off  below  300  Hz.  therefore  no  aliasing  should  have  occurred.  These 
scaled  versions  of  the  first  and  second  moments  of  the  received  optical  power  S  are  transmitted  across  an 
IEEE-488  bus  to  the  IBM-AT  computer  controller.  A  Hitachi  6301  microprocessor  controls  the  interface 
and  the  triggering  of  requests  from  the  controller.  The  digital  processing  could  have  been  done  by  the 
IBM-AT  alone:  however,  the  AT  was  intended  to  control  another  dev  ice  simultaneously .  and  it  was  therefore 
determined  to  offload  the  easy,  but  software  time-consuming  task  to  a  dedicated  board.  The  input  to  the 
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A/D  (the  low-pit-  tiltered.  scaled  version  of  the  PMT  output  I  was  also  sent  loan  oscilloscope,  which  provided 
a  useful  diagnostic  aid.  B;.  triggering  the  electronics  unit,  the  AT  returned  the  sum  of  the  samples  (denoted 
-x  *  aid  the  sum  squared  of  the  samples  (denoted  ix\ 

Based  on  N  =  1000  samples,  the  AT  first  calculates  the  normalized  variance  defined  as 


:i>x 


0/(0,)  = 


>  U1  1UIIVV 


mean" 


i^i 


1  y- « 

5loX'-b 


1  =  1 


(5  i 


In  t hi v  equation,  a  is  the  A/D  scaling  constant,  and  h  is  the  skv  background  level  which  must  be  subtracted 
from  the  mean  since  it  is  unrelated  to  the  starlignt.  and  otherwise  would  decrease  the  \ anance.  The  v  alue  oi 
b  is  computed  as  the  arithmetic  mean  ot  five  samples  of  the  background  lev  el  taken  iust  out  of  the  fieid-ot- 
v  tew  ot  the  star  to  be  tracked  Next,  the  raw  isoplanatic  angle  H  i  cf>  .A.  i  is  computed  b\  scaling  the  normalized 


STAR 


\  v 
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variance.  Finally .  the  raw  isoplanatic  angle  is  converted  to  a  zenith-corrected  isoplanatie  angle  0  (O./u  by 
dividmg  the  raw  isoplanatie  angle  by  sc.  '  (<h  ). 

IV.  VERIFICATION 

The  isoplanometer  \\  as  calibrated  in  tw  o  stages  —  first  electrically  and  then  optical  ’y.  First,  an  HP7562  A 
dynamic  signal  analyzer  w  as  used  to  electrically  verify  the  transimpedance  amplifier/voltage  amplifier/filter 
combination  in  the  electronics  unit.  In  order  to  emulate  the  PMT  output,  a  i  -MO  series  resistor  was  used  to 
convert  the  signal  generator  voltage  source  to  a  current  source.  Gams  and  offsets  were  trimmed  as 
appropriate  The  current-to-voltage  transfer  function  of  the  electronics  unit  is  shown  in  Figure  j.  Tic  7-dB 
point  is  sufficiently  high  to  pass  near’y  all  stellar  intensity  fluctuations.  Th.s  low -pass  filter  limited  PMT 
noise,  and  reduced  sample  aliasing.  Noise  measurements  verify  tha’  PMT  shot  noise  is  well  below  sampling 
resolution. 

Nexi.  a  function  generator  was  connected  to  the  electronics  unit,  and  a  square  wave  of  know  n  amplitude, 
frequency  ,  and  DC  offset  was  input.  It  was  verified  that  tne  isoplanometer  sy  stem  properly  computed  the 
mean  variance,  and  normalized  vai  mce  statistics  of  this  waveform. 

Finally .  the  setup  shown  in  Figure  6  was  used  to  perform  an  equivalent  opcca!  calibration.  A  HeNe  User 
was  sent  through  a  50/50  chopper  wheel,  and  the  chopped  beam  was  directed  through  the  telescope  and  on 
to  the  PMT.  The  A/D  input  signal  from  the  ele  ironies  unit  was  observed  on  the  oscilloscope,  and  it  was 
verified  that  the  correct  statistics  were  computed 

V.  THE  EXPERIMENT 

Unfortunately,  the  instrument  tests  originally  planned  for  Mt.  Wilson  were  canceled  due  to  the  LITE 
program  funding  cuts.  A  comparison  run  of  our  instrument  w  ith  Professor  Don  Walters'  isoplanometer  was 
similarly  canceled."'  Uncooperative  weather  conditions  in  Uexington.  Massachusetts  during  the  month  ol 
July  permitted  only  three  evenings  of  measurements,  two  of  which  were  required  to  debug  the  -etup.  Data 
collected  on  the  remaining  evening  of  25  July  IdSS  are  reported  and  analyzed  in  Section  VI. 

The  experimental  procedure  was  as  follows.  First,  the  telescope  was  set  up  on  the  A-building  loot  at 
Uincoln  Uaooratory .  lev  eled,  and  polar  aligned.  Polar  alignment  was  close  enough  so  that  the  Celestron  right- 
ascen  -ion  tracking  drive  maintained  track  over  the  duration  of  each  measurement.  Starting  with  a  bright  star 
(usually  Arcturus  or  Vegat.  the  PMT  bias  voltage  was  set  at  approximately  700  V.  For  each  star  measured, 
.he  actual  voltage  was  adjusted  so  that  the  mput  to  the  A/D  averaged  about  4  V.  as  observed  on  the 
oscilloscope.  The  oscilloscope  also  verified  that  the  peak  intensity  fluctuations  did  not  exceed  the  saturation 
voltage  of  the  A.D  <10  V).  The  IB.M-AT  software  was  started,  and  the  relevant  parameters  were  input 
manually .  The  program  then  requested  that  the  telescope  field-of-view  be  set  to  the  background,  where  five 
background  samples  were  taken,  averaged,  and  stored.  The  telescope  was  repositioned  on  the  star,  and  the 
measurement  process  began.  The  1BM-AT  recorded  local  time,  background  level,  mean  ai  d  variance  from 
the  electronics  unit,  raw  and  corrected  isoplantic  angles,  and  the  zenith  angle. 

Nominally.  5  min  of  continuous  data  were  taken  per  star  measurement.  Each  trial  consisted  of  N  =  1  OOP 
samples  taken  at  a  1  -kHz  rate.  U.r  a  total  of  1  s  oi  data.  The  system  could  manage  approximately  50  trials 
per  minute.  A  subtle  problem  in  the  electronics  unit  prevented  continuous  measurement.  Once  a  set  of  trials 
(usually  200  or  250)  was  completed,  another  star  was  chosen  and  the  process  was  repealed.  We  switched 
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between  a  single  pair  of  stars,  having  different  zenith  angles,  approximately  once  every  10  min  Our  goal 
was  to  verify  tr.e  zenith-angle  dependency  described  in  Section  II  [Equation  (4)].  On  the  evening  of  25  July 
14SS.  it  was  possible  to  complete  approximately  50  runs  between  9:30  p.m.  and  midnight. 


VI.  RESULTS 


The  data  taken  on  the  evening  of  25  July  198K  are  given  in  Figures  7  througn  1 1.  Measurements  were 
alternated  between  the  stars  Arcturus  and  Vega,  switching  approximately  every  10  min.  From  9:30  p.m  to 
1 2:00  p.m.  local  time,  the  zenith  angle  for  Vega  went  from  241'  to  1  1".  w  ith  a  minimum  of  0"  occurring  about 
1 1 .00  p.m.  Arcturus  went  from  37''  to  70"  during  the  same  period. 

Figure  7  show  s  tiie  average  zenith-corrected  isoplanatic  angle  for  each  10-min  measurement  session,  li 
local  atmospheric  conditions  are  constant,  then  the  average  values  on  consecutive  measurements  should  be 
identical  for  the  two  stars.  However,  this  j*.  not  the  case  in  Figure  7.  except  early  in  the  evening  when  both 
stars  have  small  zenith  angles  (less  than  30" i.  There  are  two  possible  explanations  for  this  discrepancy .  One 
is  that  the  turbulence  on  that  evening  might  not  nave  been  direction  independent  (i.e..  C‘  was  a  function  ol 
location  and  height,  rather  than  just  location  i. 

A  more  plausible  explanation  is  that  the  turbulence  was  strong  and.  consequently,  the  normalized 
variance  was  saturated  In  this  situation,  the  isoplanometer  theory  presented  in  Section  II  is  no  longer  valid, 
and  Equo’ion  (4)  cannot  be  used.  Note  that,  for  large  zenith  angles  in  strong  turbulence,  the  normalized 
variance  will  saturate  and.  thus,  the  isoplanatic  angle  will  be  overestimated  (as  in  Figure  7i.  Previous 
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isoplanometers  built  by  Walters  begin  to  saturate  at  normalized  variances  of  about  0.05  (see  p.  86  in 
Reference  3 1.  While  tracking  the  star  Vega  normalized  variances  were  about  0.02  to  0.03.  However,  w  hile 
tracking  Arcturus.  normalized  variances  were  0.06  or  above  after  the  initial  data  set.  Note  that  the  isoplanatic 
\  alues  for  Arcturus  are  higher  than  those  for  Vega:  under  saturated  variance  conditions,  the  isoplanatic  angle 
w  ill  be  overestimated.  It  is  plausible  that  the  instrument  was  correctly  calculating  isoplanatic  angle  early  in 
the  evening,  and  later  for  Vega  only.  Arcturus  values  later  in  the  evening  were  saturated,  and  therefore 
incorrect.  However,  this  cannot  be  stated  as  a  certainty . 

Figures  8  and  9  show  the  complete  isoplanatic  angle  data  collected  for  Vega  on  the  evening  of  25  July 
1988.  Figures  10  and  1 1  show  similar  plots  for  Arcturus  on  ths  same  night.  Note  that,  on  these  plots, 
"dropouts"  are  apparent.  These  represent  tracking  errors  t  telescope  not  stay  ing  aligned  on  the  star )  or.  in  some 
cases,  clouds  interfering  w  ith  the  measurement.  Such  effects  were  taken  out  by  post-processing,  and  arc  not 
reflected  in  the  analysis  here. 

Figures  12  and  13  plot  the  normalized  variances  vs  the  secant  of  the  zenith  angle  on  a  log-log  scale,  as 
suggested  by  Walters.10  for  both  data  sets,  and  for  Arcturus  only.  The  slope  of  the  linear  fit  is  below  the 
predicted  8/3  value,  a  result  consistent  with  a  saturated  normalized  variance. 

Finally .  Figure  14  show  s  the  saturation  isoplanatic  angle  as  a  function  of  zenith  angle.  Data  points  abo\  e 
or  to  the  left  of  the  saturation  line  are  considered  valid;  those  below  or  to  the  right  of  the  line  are  considered 
not  valid.  Note  that  all  data  for  Arcturus  are  within  the  invalid  region,  while  data  for  Vega  are  either 
unsaturated  or  marginal.  Therefore,  comparing  estimates  of  isoplanatic  angle  made  with  Vega  with  those 
made  from  Arcturus  is  invalid.  Also  note  that  if  a  fully  operational  instrument  existed,  this  technique  could 
be  used  to  predict  poor  atmospheric  conditions  l  i.e..  as  low  as4-p.rad  angle  at  0.5  pm ).  as  long  as  a  reasonably 
bright  star  existed  at  zenith  angle  less  than  38". 

VII.  CONCLUSIONS 

Group  67  has  designed  and  built  an  isoplanometer.  The  theory  of  this  device  is  based  on  the  well-known 
scintillation  averaging  approach.  The  device  was  operated  in  Lexington.  Massachusetts  for  verification 
purposes.  Lexington  is  a  poor  astronomical  site,  and  often  exhibits  a  high  level  of  turbulence.  Consequently . 
isoplanatic  angles  can  be  very  low  and  the  normalized  variance  may  saturate  under  almost  all  conditions  in 
Lexington.  The  instrument  was  designed  to  accurately  measure  isoplanatic  angle  when  the  angle  is  aboxe 
approximately  5  to  6  p.rad  at  5000  A.  the  approximate  requirement  for  the  LITE  experiment.  Initial  plans 
for  device  verification  included  co-location  with  another  isoplanometer  at  a  high-quality  astronomical  site, 
where  isoplanatic  angle  levels  would  reach  7  to  10  p.rad  at  5000  A.  The  wavelength  scaling  to  0.86  pm  (the 
LITE  operating  wavelength)  would  be  a  factor  of  1.919.  giving  isoplanatic  angles  above  14  p.rad  at  this 
wavelength  at  a  good  site.  Unfortunately,  this  extensive  verification  was  not  possible  and.  therefore, 
verification  of  dev  ice  operation  was  not  completed.  However,  results  indicate  expected  operation  under  the 
given  conditions  in  Lexington.  Verification  testing  of  this  instrument  would  therefore  require  that  it  be  placed 
at  a  "good"  astronomical  site,  such  as  the  Mount  Wilson  Observatory. 

It  should  be  noted  that  Rein  Teoste  and  Jim  Daley  of  Lincoln  Laboratory  's  Division  5  have  performed 
isoplanatic  angle  measurements  at  Firepond  using  an  instrument  built  by  Don  Walters.  They  reported" 
isoplanatic  angles  on  the  order  of  5  pLrad. 
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